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Abstract 

We introduce a model for a two configurations system, and we study the 

transition from quantum to classical behaviour. We first consider the effect of 

the interaction with the environment as an external noise and we show that 

it produces decoherence and suppression of tunnelling. These features are 

widely accepted as definition of classicality, while we believe that classicality 

implies that quantum delocalized states spontaneously evolve into localized 

ones. We than show that this evolution take place only when both noise and 

non linearity in the equations are present. 
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In this paper we focus our attention on a very old question: how and why macroscopic 
objects behave classically. If we disregard any answer which invokes a non-quantum or 
quantum-modified status for macroscopic objects we remain with the widely accepted 
point of view that environment induces decoherence. This fact, which has been shown to 
hold in a huge number of models (see for example ]§-[7|), can be stated as follows: the 
density matrix corresponding to a superposition of eigen-states of a macroscopic variable 
becomes very soon almost diagonal. Therefore, the superposition is for all practical purposes 
indistinguishable from the statistical mixture, whose probabilities are given by the diagonal 
elements |§. 

Our point is that this is not a complete answer, since it has only a statistical meaning, it 
remains to be explained why individual macroscopic objects are always found in a localized 
state. If we want such an explanation we are left with only a possibility: to show that 
any superposition quantum state spontaneously evolves in a eigen-state of the macroscopic 
variable. 

We study the problem by considering a system with two minimal energy configuration 
positions. Systems of this kind are, at low temperature, classically localized in one of the 
two minima. On the contrary, quantum behaviour allows for delocalized states and coherent 
tunnelling. Real objects of this type are, for example, superconducting rings with a weak 
junction crossed by a magnetic flux or a class of tetraedrical molecules. 

We start considering ah initio a two configurations model. The isolated system tunnels 
coherently and periodically. Then we consider the effect of a time dependent perturba- 
tion (interaction with the environment). Now there is decoherence and, for strong noise, 
suppression of tunnelling (see also |9|-|Tl|). 

Our point, as explained before, is that decoherence and suppression of tunnelling which 
implies stability of localization are not enough for classicality. Classicality needs spontaneous 
evolution into one of the two possible localized states. 

Therefore, in a next step we consider the same model, but we assume a non linear 
evolution. In absence of perturbation, it turns out that the system behaves qualitatively 
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as in the linear unperturbed case. But, this time, if interaction with the environment 
is introduced, one has that there is spontaneous evolution into one of the two localized 
states. The conclusion is that both non linearity and noise are fundamental ingredients for 
classicality. 

The density matrix can be written in a convenient form using the vector x = (x, y, z) 
defined by 



P= 2 



( 1 • \ 

1 + z x — xy 

(1) 

^x + iy I — z ^ 

The trace of this density operator equals 1 (normalization) and the determinant vanishes for 
a pure state. This last requirement implies |x| = 1 (a pure state is represented by a point 
on the unitary sphere). 

The third component z of this vector encodes the informations about the localization 
(z = ±1 means that the system is exactly localized in one of the two states) while x and y 
encode the informations about coherence (when they vanish there is complete decoherence) . 

The above considerations show that the quantum evolution p = —i[H,p] corresponds to 
the motion of a point on a unitary sphere. When the point is in one of the two poles the 
system is localized. 

The hamiltonian which describes an isolated two configurations system allowed for tun- 
nelling is typically 

H = aa x (2) 

where a x is the Pauli matrix. This hamiltonian produces coherent quantum tunnelling of 
period ^ between the two macroscopic configurations of the system (in mathematical terms, 
the two eigen-states of the Pauli matrix a z ). 
The corresponding equations for x are 

x = 

y = —2az (3) 
z = 2ay 



It is clear that (0) induces a rotation with constant angular velocity around the x axis. 
The value of the first component of the vector remains unchanged. Therefore, (0) never 
allows for localization or decoherence. 

Assume now that the interaction with the environment is described by 

H = aa x + e{t)a z (4) 

where a x and a z are Pauli matrices and e(t) is a a given realization of a stochastic process. 
The presence of a non vanishing e(t) means that we introduce an interaction with the 
environment which randomly breaks in time the energy symmetry of the two configurations. 
We find that the vector x satisfies 

x = -2e(t)y 

y = -2e(t)x - 2az (5) 
z = +2ay 

which is the superposition of two rotations, the first around the x axis with constant angular 
velocity 2a, the second around the z axis with time dependent angular velocity 2e(t). 

Assume e(t) = (3i](t) where rj(t) is a given realization of a white noise (i.e. w(t) = 
jQT](s)ds is a brownian motion). In Stratanovitch notation the only thing we have to do is 
to substitute e(t) with (3rj{t) in the above equations. Nevertheless, it is much more practical 
to use Ito notation for which the above equations (||) rewrite as 

dx = —2/3 2 xdt — 2f3ydw 

dy = -2f3 2 ydt - 2azdt + 2(3xdw (6) 
dz = +2aydt 

From (^|), taking the expectation value, we immediately obtain 

i = -2(3 2 x 

V = -2(3 2 y - 2az (7) 
z = +2ay 
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Since the density matrix depends linearly on x, the above equation describes the evolu- 
tion of the averaged density matrix p. This matrix describes the mixture of all the states 
associated to the ensemble of all the realizations of the noise in the hamiltonian (|J). 

Equation (|7|) can be easily solved. When (3 2 < 2\a\ the solution is: 



x (t) =e- 2 r*x(0) 

^(t) = e -P 2t [y(0) cos(wt) + ci sm(ut)] (8) 



z{t) = e~P *[z(0) cos(W) + c 2 siii(wt)] 
where u = \/\P 4 — 4a 2 1, c\ = -/3 y ^~ 2az ^ an d c 2 = - z ^ +2ay ^ . Looking at the third com- 



ponent z(t) one realizes that in this region a quantum coherent behavior partially survives 
to the noise for a certain time. 

When (3 2 > 2\a\ the solution is purely exponential and is given by 



x(t) =e" 2/3i x(0) 

^(t) = e - p2t [y(0) cosh(wt) + ci sinh(o;t)] (9) 



z{t) = e'P '[^(0) coshfW) + c 2 mnh(ti^)] 
In this strong noise region the coherent behaviour is completely destroyed since the 



localization probability relaxes exponentially. The interesting fact is that for large /3, z(t) 
relaxes with an exponent which vanishes as This fact means that the time necessary for 
derealization becomes infinite in the limit of extremely strong noise. 

Summarizing, we have damped oscillations for a weak interaction with the environment 
(f3 2 < 2 1 a |) and incoherent relaxation for strong interaction ((3 2 > 2\a\). Furthermore, in 
the limit of very strong noise a localized configuration turns out to be almost stable. 

Since in both cases x — >• exponentially, we have that 

which says that in any case we have almost complete decoherence (no diagonal terms) at 
large times, the non-diagonal elements being damped whereas the diagonal ones carry the 
probabilities. 



From (|T0|) we could conclude that an individual system is in one of the two configuration 
positions. This is not true, in fact, from the stochastic equations (0) one can derive the set 
of linear equations 

z 2 = 4azy 

f = -Aazy - Afy 2 + 4{3 2 ^ (11) 
yz = —2/3 2 zy + 2ay 2 — 2az 2 

Using the condition x 2 + y 2 + z 2 = 1, one can replace x 2 with 1 — y 2 — z 2 in the second of 
the above equations. For any j3 ^ one finds that the solution converges to the stationary 
solution x 2 = y 2 = z 2 = |, z~y = 0. This means that the typical state is not localized 
(localization for all states would imply z 2 = 1). A new ingredient is thus necessary to 
produce localization. 

The effect of the interaction with the environment can be taken into account also con- 
sidering a non linear term in the differential equations that describe the evolution of the 
density matrix. In our case, we can assume that a depends on the state (a = a(x,y,z)). 
For the sake of simplicity assume, as usual, that this dependence reduces to a dependence 
on the modulus square of the wave function. In our language it means that a depends only 
on z. There are many reasons for this choice, that we will discuss in a forthcoming paper 
0- 

Also assume that: 

a) a(z) is an even function of z: a(z) = a(—z); 

b) a(z) is zero at the poles: ct(+l) = a(— 1) = 0, positive otherwise and sufficiently 
smooth around the poles. 

The requirement b) implies that a localized state remains, in fact, localized. Nevertheless, 
in absence of noise, the system is not able to spontaneously localize. In this case, in fact, 
equation (|3|) still holds, with the only difference that a depends on z. 

Notice that the first component x of the vector x remains constant during the motion. 
This means that there is no decoherence, and, furthermore, that the system can never reach 

6 



one of the poles. Indeed, the motion remains periodic, since we are on a two dimensional 
variety. The only exception is when the system is initially on the meridian x = 0, in which 
case it moves along this meridian toward one of the poles, producing localization, but this 
ensemble of initial conditions has probability zero. 

In conclusion, non linear differential equations, like the equation modelizing the inter- 
action with the environment, are not able to reproduce a classical behaviour in the more 
restricted sense of induced localization. Nevertheless, this goal can be achieved if one takes 
into account the simultaneous presence of these two different mechanisms. 

Let us consider again the hamiltonian (|4]) with a = a(z); the previous requirements on 
a(z) are also assumed. The differential equations (|6|), which remain unchanged, have two 
fixed points corresponding to the poles. What we can show is that, independently on the 
initial conditions, one has z(t) — > ±1. This convergence is not due to the attractiveness 
of the poles but to the fact that the motion along the meridians becomes more and more 
slow approaching the poles. In other words, the system stays for long time intervals around 
them, the distribution of this time intervals having sufficiently long tails to guarantee that, 
at the end the system will spend almost all its time around a pole (it will be almost surely 
in a pole). 

We can sketch how to show that the system localizes for large times; a rigorous proof 
being beyond the scope of this paper. The requirements on a(z) imply that it can be 
approximated around the poles by 



We now prove that the limiting (steady) distribution is concentrated on the poles. Using 




(12) 



equations @ and taking the average, we have 





(13) 
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where < z c < 1, and #(•) is the step function. 

We can safely assume that any initial distribution on the surface of the unitary sphere 
will evolve, after a transient time, toward a stable distribution. For this stable distribution, 



the above time derivative of -^0(z — z c )y vanishes. 
Furthermore, from equations (|6|) it also follows that 

d 



dt 



6(z-z c ) = 2a(z)5(z - z c )y (14) 



The above equation implies for the steady distribution 5(z — z c )y = 0, which, in turn, implies 



6(z - z c )y = 0. 

This last equality, inserted in ( pT3|) gives for the stationary distribution 



a(z)6(z - z c )z = a(z c )5(z - z c )y 2 (15) 



which leads to 



lim a(z)9{z ~ Zc)z = 1 (16) 
z ^ a{z c )5{z - z c )y 2 

It is intuitive, and it can be also easily shown (again from (^|)) that the relative difference 



between S(z — z c )y 2 and S(z — z c )x 2 vanishes in the limit z — > z c , in fact the random rotation 
around the z axis becomes very fast compared with the rotation around the x axis. Therefore, 



since x + y + z = 1, we can replace S(z — z c )y 2 with S(z — z c )^- in equation (|i6|), which 
we rewrite as 

j!(l-z 2 )zp(z)dz 



lim JZc \ 1 ' rv ' = 1 (17) 

where we have used the expansion of a(z) ~ chq(1 — z 2 ) around the north pole and we have 
introduced explicitly the steady state probability p(z) that the system is at a quote z. 

Assuming that at p(z) ~ ^ around the north pole with 7 < 1, we find that the 
above equation gives 

2^ = 1 < 18 > 



which cannot be satisfied for any normalizable distribution. 

Finally, since ( |T5| ) is satisfied for a distribution which is concentrated in the north pole 
we find that this is the only possible stable distribution. Repeating the argument for — 1 < 
z c < 0, and taking into account the symmetry of the system (|6]) with respect to the changes 
y — > —y and z — > —z, one can conclude that the steady distribution is made of two Dirac's 
delta distributions centered on the poles, with equal weights. 

Notice that this result implies again (|l0f) . Therefore, for what concerns decoherence 
nothing is changed, the important difference being that now p 2 = ~p for large times. 

In fig. 1, where z(t) is plotted starting from a numerical solution of the differential 
equations @ for a = 1 and (3 = 7, one sees clearly that the system spends a large part of 
the time in very narrow regions around the poles. 

The time necessary to make a transition between the poles turns out to be negligible with 
respect to the long periods in which z ~ ±1, and therefore, if one perform an observation 
of the system, it is almost sure to find it in a localized configuration. A useful quantity in 
order to show this fact is the time average of z(t) 2 , defined by: 

L(t) = 1 f[z{s)fds (19) 
t Jo 

This average tends to grow in time since the periods when z(t) 2 is substantially different 
from 1 become soon negligible, as shown by fig. 2, where L(t) is plotted for the same noise 
realization as fig. 1. 

In conclusion, we have presented a model where the coupling with the environment is 
represented both by a noisy potential and a non linear correction to the differential equations. 
This second effect should be a consequence of the feed-back of the system which is able to 
encode informations on the macroscopic system. Non linearity and noise are not able to 
produce localization, while together they induce the system to localize in one of its minimal 
energy configuration positions. Moreover, the localization positions are chosen randomly 
and they are not fixed in time, in the sense that the system can always make a transition 
between the two different macroscopic states. Therefore, it is not possible to predict the 
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result of an observation. 

As a final remark, we would like to stress that our model is based on a genuine quantum 
dynamics without more or less phenomenological dissipative terms in the hamiltonian [|I],0 . 

In our opinion the results presented in this paper could provide a general framework for 
describing the subtle transition associated with the emergence of the classical world. 

It is a pleasure for us to thank Gianfausto Dell'Antonio for fruitful discussions and useful 
suggestions. 
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FIG. 1. A typical realization of z(t) as function of t for a® = 1 and (3 = 7. The differential 
equation [6] is numerically solved with a time step of 10 -2 . 
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FIG. 2. L(t) (19) as function of t for the same noise realization of fig. 1, with ao = 1 and (3 = 7. 
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